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1. Introduction 

As control of quantum phenomena has become an increasingly realistic objective 
due to both theoretical and technological advances [Q, technical issues such as 
the controllability of quantum systems have received considerable attention recently 

i I, i 1, 1, 0, 1, 1, |io|,0,|il,|il. 



This has led to the introduction of various definitions of controllability for quantum 
systems, such as complete controllability, wavefunction controllability, density matrix 
controllability, observable controllability, etc., as well as a search for necessary and 
sufficient conditions to characterize these different degrees of controllability. Using 
general results from Lie group and Lie algebra theory, necessary and sufficient criteria 
for various types of controllability have been derived ||10|, |14| and for certain types 



of quantum control systems such as A^-level ladder systems with dipole interaction 
1^, ^ or interacting spin systems [Tl|, |15[, it has actually been demonstrated that 



these criteria are satisfied. Another interesting result is the graph connectivity theorem 



established in |T2[, which provides sufficient conditions for wavefunction controllability 
for A^-level systems, if the energy levels are non-degenerate and the transition frequencies 
are distinct, as well as generalized versions applicable to other notions of controllability 



In this paper, we consider the implications of these theoretical results for atomic 
systems with degenerate energy levels and transition frequencies. The remainder of 
the paper is organized as follows. In section two, we review basic concepts and we 
give precise definitions for various degrees of controllability. Section three provides a 
concise summary of necessary and sufficient conditions for each type of controllability. 
In section four these results are applied to determine the degree of controllability for 
atomic transitions between two energy levels with varying degrees of degeneracy. The 
implications of the results, as well as limits of this approach and directions for future 
work, are discussed in sections five and six, respectively. 

2. Basic concepts 

2.1. Representation of Quantum States 

The state of any A^-level quantum system can be represented by a positive-definite, 
trace-one operator p acting on a Hilbert space. The representation of the state is unique 
under the conditions described in [Q, chap. 5. This so-called density matrix always 



has a discrete spectrum with non-negative eigenvalues Wn that sum to one, Y.n'^n = 1, 
and a spectral resolution 

N 

p=Y,Wn\^n){'^nl (l) 
n=l 

where |\Ef„) are the eigenstates of p The eigenstates corresponding to different 



eigenvalues are always orthogonal and even if an eigenvalue occurs with multiplicity 
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greater than one, the can always be chosen such that they form an orthonormal 
set in the Hilbert space Ti. Moreover, by adding states with eigenvalue w„ = if 
necessary, we can always assume that {|^n), 1 < < forms a complete orthonormal 
set, i.e., a basis for the Hilbert space Ti. The are the corresponding dual states 
defined by 

I = 5mn Vm, n. (2) 

Definition 1 A density matrix p represents a pure quantum state if it has rank one. 
In any other case p represents a mixed quantum state. 

If p represents a pure state then its eigenvalues counted with multiplicity are 
{1,0, ...,0}, i.e., there is one non-zero eigenvalue Wi = 1 and all other eigenvalues 
are zero. Mixed states must not be confused with superposition states. For a physical 
system there usually exists a preferred basis consisting of the energy eigenstates of 
its internal Hamiltonian. Superposition states are simply linear combinations of these 
energy eigenstates. However, they are pure states, as are the energy eigenstates. 

Example 1 Consider the density matrices 

/ i \ 





V i / 

Pi represents the basis state while p2 represents the superposition state |^) = 
;^(|1) + |4)) = -^{1,0,0,1)'^ since we have p2 = |\E')(\E'|. However, both density matrices 
have eigenvalues {1,0,0,0} and thus represent pure states, pz, on the other hand, 
represents a true mixed state since we have ps — ||1)(1| + ||4)(4|, i.e., its eigenvalues 
are {|, |, 0, 0} and ps has rank two. 



Pi 



/ 1 \ 





V / 



P2 



/i i\ 





U 2 



P3 



2.2. Kinematical Equivalence Classes of States 

Conservation laws such as the conservation of energy and probability require the time 
evolution of any (closed) quantum system to be unitary. Thus, given a pure state |^o)) 
its evolution is determined by — ^(Ol^o) where U{t) a unitary operator for 

all t and f/(0) = /. Consequently, a general quantum state represented by a density 
operator po must evolve according to p{t) = U{t)f)oU{ty with U{t) unitary for all times. 
This constraint of unitary evolution induces kinematical restrictions on the set of target 
states that are physically admissible from any given initial state. 

Definition 2 Two quantum states represented by density matrices po and pi are 
kinematically equivalent if there exists a unitary operator If such that 

pi = UpoUl (3) 
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Theorem 1 Two density matrices po and pi are kinematically equivalent if and only if 
they have the same eigenvalues. 

Proof: Given two N x N density matrices that have the same set of eigenvalues, 
we can always find two sets of orthonormal eigenvectors and a unitary operator 
U G U{N) that maps one set of eigenvectors onto the other such that |\E'n) = f/|0„) for 
all n. Thus, we have 

Pi = ^W„|^„)(^„| = ^Wnt>|0n)(0n|f/^ = U PoU^ ■ 

n n 

Similarly, if holds then pi and p2 must have the same eigenvalues since suppose 

Po = En«^n|0n)(0n| then 

Pi = UY.^n\<l>n){<Pn\U^ = \<l>n) {<Pn\U^ 

n n 

Since U is unitary, f/|0n) is an orthonormal basis for pi and the eigenvalues of pi are 



2.3. Dynamical Lie Groups and Reachable States 

For any given initial state, only states in the same kinematical equivalence class can 
possibly be dynamically reachable. However, the set of dynamically reachable states 
may be a subset of the kinematical equivalence class. In general, the set of states that 
are dynamically accessible from a given initial state depends on the dynamical Lie group 
of the system. Consider a quantum system with a control-linear Hamiltonian 

M 

H[f,it), . . . , /M(t)] = i^o + E Ut)Hm, (4) 

m=l 

where the fm, ^ ^ ^ M, are (independent) bounded measurable control functions. 
Since the time-evolution operator U (t) has to satisfy the Schrodinger equation 

thjUit) = H[f,{t),...jMit)]Uit), (5) 

only unitary operators of the form 

U{t) = exp+ l^-^H . . . , /M(t)]} , (6) 

where exp_,_ denotes the time-ordered exponential, qualify as evolution operators. Using 
the Magnus expansion of the time-ordered exponential, for instance, it can be seen that 
only unitary operators of the form exp(a;), where x is an element in the dynamical Lie 
algebra C generated by the skew-Hermitian operators iHo, ■ ■ ■ j^Hm, are dynamically 
realizable. These unitary operators form the dynamical Lie group S of the system. 

For practical purposes, it is often more convenient to consider the related Lie algebra 
C, which is generated by the trace-zero parts of the operators iHm, 

fr _fr MHm) f 

J^m — J^m j;^ J-Nt \l) 
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for < m < M, where 1^ is the identity matrix of dimension A^. C is always 
a subalgebra of su{N) since it is generated by trace-zero, skew- Hermit ian matrices. 
Furthermore, if Ti{Hm) = for all m then L = L\ otherwise we have C = L ® u{l) 
where u{l) is spanned by the identity matrix Jtv. 



3. Degrees of controllability 

3.1. Definitions 

In this section we give precise definitions for various degrees of controllability for 
quantum systems, which are relevant in areas such as quantum-state engineering. 



quantum chemistry |12| and quantum computing |18 



Definition 3 A quantum system is completely controllable if any unitary evolution 
operator U is dynamically accessible from the identity I . 

That is, a quantum system is completely controllable if there exists T > 0, a 
set of admissible control functions . . . , /^(t)) defined for < t < T and a 

corresponding trajectory U{t) satisfying the Schrodinger equation as well as the 
boundary conditions f/(0) = / and U{T) = U. 

Definition 4 A quantum system is density matrix controllable if for any given initial 
state represented by a density matrix po? '^^^ kinematically equivalent states can be 
dynamically reached at some later time T > 0. 

More precisely, a quantum control system is density-matrix controllable if given 
any two kinematically equivalent states po and pi, there exists T > 0, a set of 
admissible control functions • • • , fuif)) defined for < t < T and a corresponding 

evolution operator U (t) satisfying the Schrodinger equation (H) such that U{0) = I and 
pi = f/(T)pof/(T)t. 

Definition 5 A quantum system is pure-state controllable if for any given pure initial 
state represented by a wavefunction |\E'o), any other pure state |\E'i) can be dynamically 
reached at some later time T. 

More precisely, a quantum system is pure-state controllable if given any two 
pure states I^E'o); l^i) there exists T > 0, a set of admissible control functions 
(/i(t), . . . , defined for < t < T and a corresponding evolution operator U(t) 

satisfying the Schrodinger equation such that f/(0)|\E'o) = I^E'o) and f/(T)|\E'o) = l^i)- 

Definition 6 A quantum system is observable controllable if any observable represented 
by a Hermitian operator A onTi can dynamically assume any kinematically admissible 
expectation value for any given initial state of the system. 
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The kinematically admissible expectation values for any observable A depend on 
the initial state of the system. Precisely, the expectation value (ensemble average) of 
the observable A is bounded by |]19| 

N N 

J2 WnXw+l-n < {Ait)) < ^ WnXn, (8) 
n=l n=l 

where Wn are the eigenvalues of the equivalence class of density operators selected by 
the initial state, which are counted with multiplicity and ordered in a non-increasing 
sequence, and A„ are the eigenvalues of the operator A, also counted with multiplicity 
and ordered in a non-increasing sequence. The upper bound is assumed if the system is 
in state 

N 

p+ = (9) 

n=l 

where l^^n) are the eigenstates of A satisfying A|\l'„) = A„|\E'„). Similarly, the lower 
bound is assumed if the system is in state 

N 

p_ = J2wN+l-n\^n){^n\- (10) 
n=l 

If either A or po has eigenvalues occurring with multiplicity greater than one then 
there exists a subspace of states for which the upper and lower bounds are achieved. 
Otherwise, there is a unique state for which these bounds are realized. Any intermediate 
value can be achieved for some state pi, which is kinematically admissible from the initial 
state po- 



3.2. Necessary and sufficient conditions 

A necessary and sufficient condition for complete controllability is that the dynamical 
Lie group S of the system be U{N). Noting that S is always a subgroup of U{N), we 
see immediately that complete controllability is the strongest possible requirement of 
controllability. Since Eq. (^ defines a right-invariant control system on the compact 
Lie group U{N), or SU{N) if all the Hamiltonians Hm have zero trace, it follows from 
theorem 7. 1 in ||2^ that a necessary and sufficient condition for the dynamical Lie group 
S to be U{N) is that the dynamical Lie algebra C be isomorphic to u{N). Similarly, 
S ~ SU{N) if and only if £ ~ su{N). 

As pure states can be represented by complex unit vectors, a sufficient condition for 
pure-state controllability of a quantum system of dimension is that its dynamical Lie 
group S act transitively on the unit sphere S"^^^^ in C^. It is easy to see that transitive 
action on 5*^^^^ is necessary for pure-state controllability as well. Using classical results 
on tranformation groups of spheres |21| it can be shown that the only subgroups of 
U{N) that act transitively on the unit sphere in are U{N) itself, SU{N), and if A^ 
is even, Sp^^) well as Sp^^) ^ U{1) with f/(l) = {c^'^/tv} |TD|. Thus, a quantum 
system will be pure-state controllable if and only if its dynamical Lie group S contains 
one of the Lie groups above. In (Theorem 4) it is shown that this is possible only 
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if the dynamical Lie algebra C of the system is u{N), su{N), or if N even, or 
sp{f)®u{l). 

A necessary and sufficient condition for density matrix controllability is that the 
dynamical Lie group S act transitively on all equivalence classes of density matrices. 
Clearly, 5 ~ t/ (N) is sufficient for density matrix controllability. SU (N) is also sufficient 
since for any target state pi that is reachable from a given initial state po via a unitary 
transformation U, there exists an equivalent transformation If in SU{N) such that 
pi — UpoW . Pure-state controllability is a prerequisite for density matrix controllability 
since a pure state \^) can be represented by a density matrix |^)(^|. Thus, only systems 
whose dynamical Lie group is U{N), SU{N), S'p(y) or Sp{y) x U{1) qualify for density 
matrix controllability However, we shall see that the latter two are not sufficient for 
density matrix controllability since there are kinematically admissible mixed states that 
cannot be reached from certain initial states. 

Example 2 Let N = 2i and assume the dynamical Lie algebra C generated by the 
trace-zero Hamiltonians H,n. < m < M, is Then there exists a 2£ x 2£ matrix 

J, which is unitarily equivalent to 

such that 

{iHmfj + J{iHm)^0, Vm. 

Without loss of generality suppose J — J. Let x be a 2£ x 2£ diagonal matrix with 
diagonal {—wi, —W2, ■ ■ ■ , —w^, wi, W2, ■ ■ ■ , wi) and set 

Po^Y£^2i + x. (12) 

If < Wi < W2 < ■ ■ ■ < W£ < ^ then the diagonal elements of po are distinct, non- 
negative and sum to one. Thus po is a density matrix of rank N with N = 2£ distinct 
eigenvalues. 

ix is skew-Hermitian and satisfies {ix)-^ J + J{ix) = 0. Thus, ix is an element 
in the dynamical Lie algebra C. Hence, it has to remain in L under the action of the 
dynamical Lie group S. Thus, the state 

Pi = ^^2£ + y, (13) 

where y is a 2£x2£ diagonal matrix with diagonal {—W2, —Wi, . . . , —w^, Wi, W2, ■ ■ ■ , w^), 
which is kinematically equivalent to po, is not dynamically accessible from po since iy 
does not satisfy {iy)^J + J{iy) = and is thus not in C. 

This example shows that Sp{^) or Sp{^) x U{1) does not act transitively on all 
equivalence classes of density matrices; in particular it does not act transitively on 
equivalence classes of density matrices of rank N with distinct eigenvalues. 

Density matrix controllability is clearly sufficient for observable controllability for 
it implies that any state in the same kinematical equivalence class as the initial state 
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can be reached dynamically. Therefore, any observable can dynamically assume any 
expectation value (ensemble average) that is kinematically allowed for the equivalence 
class of states selected by the initial state. The previous example also shows that density 
matrix controllability is a necessary condition for observable controllability since pi is 
itself an observable that assumes its kinematical upper bound exactly if the system is in 
state pi, which is not dynamically accessible from the initial state po. (It is important 
here that pi and po have distinct eigenvalues since this guarantees that there is only 
a single state for which pi assumes its kinematical upper bound.) 
All the above observations can be summarized as follows. 

Theorem 2 A necessary and sufficient condition for a quantum system with 
Hamiltonian ^ and dynamical Lie algebra C (C) to be 

(i) completely controllable is C ~ u{N); 
(a) density-matrix controllable is C c::^ su{N); 
(Hi) observable controllable is C c:^ su{N); 

(iv) pure-state controllable is C ^ su{N), or if N is even, C ^ sp{y)- 

This theorem is not only theoretically interestesting. It can be used in practice to 
determine the degree of controllability of a quantum system by computing the dynamical 
Lie algebra C of the system and determining whether it is isomorphic to one of the 
algebras above. For specific model systems one can often even determine the Lie algebra 
as a function of the parameters of the model and apply the previous theorem to decide 
if the system is controllable and to what degree |]^, ^. 

Notice that density matrix and observable controllability are equivalent. Complete 
controllability is theoretically slightly stronger although the difference is subtle and can 
usually be ignored in practice. To see why this is the case, observe that if Tr(^m) = 
for all m, i.e., in particular if Tr(ifo) = then there are certain unitary operators that 
are not dynamically attainable. For instance, if £ = su{2) then the unitary operator 



is not dynamically realizable unless is an integer multiple of 2tt. Hence, given 
l^o) = ^(|1) + |2)) then the state l^^i) = ^(|1) + e*'^|2)) is technically not reachable 
from |\l>o). However, the state |^i) = -^{e~''*'^'^\l) + e''^/^\2)), which differs from 
only by an absolute phase factor e"""^/^ is dynamically reachable. Therefore, the trace 



4. Controllability of atomic systems 

We shall now apply the results of the previous sections to determine the degree of 
controllability of various atomic systems. We consider an atomic system with two main 
energy levels, where the number of degenerate sublevels depends on the F-value of the 
atomic level and is given by the formula 2F + 1. We can couple magnetic sublevels 




of the operator Hq is usually not significant for practical purposes. 
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Ex 1: F=0 / F=0 Ex 3a: F=l / F=l Ex 4a: F=l / F=2 Ex 5a: F=2 / F=2 

|2) |2) |4) |6) |1) |3) |5) |7) |8) |2) |4) |6) |8) |10) 

I III "III" mil 

|1) 11) |3) |5) |2) |4) |6) 

Ex 2: F=0 / F=l Ex 3b: F=l / F=l Ex 4b: F=l / F=2 

|2) 13} |4) |2) |4) |6) |1) |3) |5) |7) |8) 

""" m ^m:" 

|1) |3) |5) |2) |4) |6) 

Ex 3c: F=l / F=l Ex 4c: F=l / F=2 

|2) |4) |6) |1) |3) |5) |7) 




|1> |3) |5) |7) |9) 
Ex 5b: F=2 / F=2 

|2) |4) |6) |8) |10) 






|1) |3) |5) |7) |9) 
Ex 5c: F=2 / F=2 

|2) |4) |6) |8) |10) 




|1) |3) |5) |2) |4) |6) 

Ex 3d: F=l / F=l Ex 4d: F=l / F=2 

|2) |4) |6) |1) |3) |5) |7) 





|1) |3) |5) |2) |4) |6) 

Figure 1. Coupling diagrams for various examples. 



|1) |3) |5) |7) |9) 

Ex 5d: F=2 / F=2 

12} |4) |6) |8) |10) 

|1) |3) |5) |7) |9) 



with the same quantum number m using a hnearly polarized field and sublevels with 
Am = ±1 using left and right circularly polarized fields, respectively. Fig. |I| shows the 
coupling diagrams for various cases. 



4.1. Transition between two F = levels 

Example 1 in figure |l| shows the trivial case of a transition between two non-degenerate 
[F = 0) levels. The internal Hamiltonian of this two-state system is 

and the interaction Hamiltonian for a linearly polarized control field is 

where c? 7^ is the dipole moment of the transition. Note that iHo and iHi generate 
u{2) if E2 7^ —El and su{2) if E2 = —Ei. Thus, the system is always density matrix and 
observable controllable (and completely controllable in the former case) using a linearly 
polarized field. 



4-2. Transitions between F = and F = 1 levels 

Example 2 shows a diagram for transitions between sublevels of two atomic levels with 
F = and F = 1, respectively. Since the F = 1 level is three- fold degenerate, the 
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internal Hamiltonian of this four-state system is given by 
( El Q \ 



\ 









E2 







E2 








E2 



(16) 



in the standard basis. Since a linearly polarized field will only drive transitions between 
states with the same magnetic quantum number, the interaction Hamiltonian for a 
linearly polarized field is 

\ 




( 


d 

V 



d 








/ 



(17) 



Similarly, if a left or right circularly polarized field is applied, it will only drive transitions 
between states whose magnetic quantum numbers differ by Am = ±1, respectively 
Thus, the interaction Hamiltonians for a left or right circularly polarized field are 



( 
d 

V 



d 






\ 


0/ 



( 



\d 



d \ 



/ 



(18) 



It is obvious from the diagram that all three polarizations are required to couple all the 
levels and noting that %Eq^ iEi, iE2 and iE2, generate the Lie algebra w(4) (or su{A) if 
E2 — —\Ei) shows that the system is always density matrix and observable controllable 
in this case. 



4.3. Transition between two F = 1 levels 

Example 3 shows four different coupling diagrams for transitions between sublevels of 
two atomic levels with F = 1. Considering the ordering of the states chosen in the 
diagram, the internal Hamiltonian is 
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(19) 
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and the interaction Hamiltonians for linearly, as well as left and right circularly polarized 
fields are 
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(20) 



(21) 



The Lie algebra generated by iHo and iHi is isomorphic to -^(2) (or su{2) if 
E2 = —El). Thus, if only a linearly polarized field is used then the system breaks 
up into three congruent, non-interacting two-level subsystems (as shown in Example 
3a) and behaves effectively like a two-state system. 

The Lie algebra generated by iHo, iHi and iH2 is isomorphic to sp{3) ® m(1) (or 
sp{3) if E2 — —El). Thus, if linearly and left circularly polarized fields are used then 
the system is pure-state controllable (see Example 3b) . The same holds if hnearly and 
right circularly polarized fields are used instead (not shown). 

However, even if linearly as well as left and right circularly polarized fields are 
used (as shown in Example 3c) the system is neither density matrix nor observable 
controllable since the Lie algebra generated by iHo, iHi, iH2 and iH^ is still isomorphic 
to sp(3) © u{l) (or sp(3) if E2 = —Ei) as in the previous case. 

Finally, if only left and right circularly polarized fields (see Example 3d) are used 
then the system decomposes into two congruent, non- interacting three- level subsystems 
spanned by the states |4), |5) and |2), |3), |6), respectively. Furthermore, the Lie 
algebra generated by iHo, iH2 and iH^ is u{3), i.e., the system behaves effectively hke 
a completely controllable three-state system. 



4.4- Transitions between F = 1 and F — 2 levels 

Example 4 shows four different coupling diagrams for transitions between sublevels of 
two atomic levels with F — 1 and F — 2, respectively. Given the ordering of states 
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chosen in the diagram, the internal Hamiltonian is 
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(22) 



and the interaction Hamiltonians for hnearly as well as left and right circularly polarized 
fields are 
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(23) 



(24) 



(25) 



If only a linearly polarized field is used (as shown in Example 4a) then the Lie algebra 
of the system is again u{2) (or su{2) if E2 — —Ei), i.e., the system behaves effectively 
like a two-state system. 
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If linearly and left circularly polarized fields are used then states |1) through |7) form 
a density matrix and observable controllable subsystem since the Lie algebra generated 
by iHo, iHi and iH2 is isomorphic to u{7) (or su{7) in the trace- zero case). 

If linearly as well as left and right circularly polarized fields are used (as shown in 
Example 4c) then the system is density matrix and observable controllable since iHo, 
iHi, iH2 and iHa generate the Lie algebra u{8) (or su{8) if 3Ei + 5E2 = 0). 

Finally, if only left and right circularly polarized fields are used then the system 
decomposes into two non-interacting subsystems consisting of states |3), |4), |7) and 
states |1), 1 2), |5), |6), |8), respectively, as shown in Example 4d. The system is therefore 
not controllable. However, it can be verified that each of the subsystems is effectively 
completely controllable. 



4.5. Transition between two F = 2 levels 



Example 5 shows four different coupling diagrams for transitions between sublevels of 
two atomic levels with F — 2. Given the ordering of states in the diagram, the internal 
Hamiltonian is 
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(26) 
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and the interaction Hamiltonians for linearly as well as left and right circularly polarized 
fields are 
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(28) 



(29) 



The Lie algebra generated by iHo and iHi only is isomorphic to u{2) (or su{2) if 
= —El). Thus, if only a linearly polarized field is used then the system behaves 
again effectively like a two-state system. 

The Lie algebra generated by iHo, iHi and iH2 is isomorphic to sp{5) © u{l) (or 
sp{5) if E2 = —El). Thus, if linearly and left circularly polarized fields are used then 
the system is pure-state controllable (see Example 5b) . The same is true if linearly and 
right circularly polarized fields are used instead (not shown). 

However, even if linearly as well as left and right circularly polarized fields are 
used (as shown in Example 5c) the system is not density matrix controllable as the Lie 
algebra generated by iHo, iHi, iH2 and iH^ is isomorphic to sp{5) ® u{l) (or sp{5) if 
E2 = —El). 

If only left and right circularly polarized fields are used then the system decomposes 
into two congruent, non-interacting five-level subsystems |4), |5), |8), |9) and |2), 
|3), |6), |7), |10) (see Example 5d). Since the Lie algebra generated by iHq, iH2 and 
iHs is isomorphic to m(5), the system behaves effectively like a completely controllable 
five-state system. 



5. Interpretation of results and discussion 



The main results of the analysis in the previous section can be summarized as follows. 
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(i) A transition between two F = k levels with = 1, 2, 3 is pure-state controllable 
using linearly and either left or right circularly polarized fields, but it is not density 
matrix or observable controllable even if linearly, left and right circularly polarized 
fields are applied. (The F = case is an exception — a transition between two 
F = levels is clearly always density matrix and observable controllable using only 
a linearly polarized field.) 

(ii) A transition between two levels with F = k and F' = A; + 1 for = 0, 1, 2, on the 
other hand, is always density matrix and observable controllable using linearly as 
well as left and right circularly polarized fields. 

(iii) None of these transitions are controllable using only left and right circularly 
polarized fields. 

(iv) If only linearly polarized fields are applied then all the above transitions behave 
effectively like two-state systems. 

5.1. Practical implications of the results 

As regards the implications of these results, the last observation is important because 
it provides a mathematical justification for the quite common practice of treating 
transitions between two energy levels like two-level systems even if the energy levels 
are degenerate. Indeed, our analysis shows that if only linearly polarized light is used 
(and the sublevel structure is not important for the application at hand) then this 
simplification is mathematically justified. 

Another interesting result is that the degree of controllability of a transition between 
energy levels with a different number of sub-levels is different from that of a transition 
between two energy levels with the same number of sublevels. In the former case, the 
degree of controllability analysis shows that every unitary operator in U{N) or SU (N) 
(where N is the total number of coupled sublevels) is dynamically realizable, while 
in the latter case the system is only pure-state controllable no matter what fields are 
applied. To better understand the implications of the latter statement, let us consider 
an example. 

Example 3 Suppose we have a system whose ground state is a P3/2 state with three 
sublevels. Assume that the initial populations of these sublevels are W-i, wq and w+i, 
where w±i ^ 0, Wq ^ and w^i + Wq + w^i = 1. Consider a transition to another 
(initially unoccupied) P3/2 (triplet) state. In this case, our analysis shows that it is 
impossible to selectively excite the population of any of the sublevels. See \Appendix Bj for 
details. For instance, we cannot selectively transfer the population of the m = sublevel 
to the m = level on the upper electronic surface without disturbing the populations of 
the other sublevels. This may not seem surprising as it would be difficult to imagine a 
pulse scheme that would accomplish such a task. However, if in the previous example 
the upper level is a state with more than three sublevels (e.g., a -D5/2 or Fj/2 state) then 
selective excitation of any one of the sublevels is possible since, in principle, we can 
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P. 



3/2 



P. 



3/2 



W^l Wq W+i 

not allowed! 



D 
P 



5/2 



3/2 



allowed! 



w+i 



Figure 2. Selective excitation of sublevel populations: not allowed for a transition 
between P3/2 states but possible for a transition between a P3/2 f^nd a -D5/2 state, 
according to degree of controllability analysis. 



Ex 3c(2): F=l / F=l Ex 5c(2): F=2 / F=2 

|2) |4) |6) |2) |4) |6) |8) |10) 




|1) |3) |5) |1) |3) |5) |7) |9) 

Figure 3. Coupling diagrams for transitions between two P3/2 and D^f2 states, 
respectively, for certain alkali atoms. The transition between m = states is prohibited 
even for linearly polarized light. Nevertheless, controllability analysis shows that the 
degree of controllability of the system remains the same. 



generate any unitary operator, including the one corresponding to selective excitation of 
any sublevel. See figure 

5.2. Effect of system modifications on controllability 

The question also arises, how small modifications of the system affect the degree of 
controllability. For instance, for alkali atoms transitions between the m = sublevels 
are prohibited for transitions between two states with the same F values, even for linearly 
polarized hght. Thus, given a transition between two P3/2 or states, for example, 
we must modify the transition diagrams in examples 3a and 5a as shown in figure ^. 
To reflect this change in the model, the values of the corresponding dipole moments in 
the Hamiltonian Hi (see equations ( pOD and (pTf), respectively) are set equal to zero. 
Nevertheless, computation of the Lie algebra shows that in both cases the Lie algebra 
remains isomorphic to sp{3) and sp{5) respectively, even with the same J matrix. Thus, 
the changes do not affect the degree of controllability of the system. 

6. Limits of controllability analysis and open problems 
6.1. Constructive control using the Lie algebra structure 

The degree of controllability of a system allows us to determine what control objectives 
are realistic, at least in principle, which is important for practical applications. As we 
have shown in the previous section, while controllability analysis based on the dynamical 
Lie algebra of the system often confirms our intuition about the controllability of the 
system, it can sometimes lead to surprising results, such as the difference in the degree of 
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controllability for transitions between states with the same F value, and those between 
states with different F values, etc. 

Precise knowledge of the dynamical Lie algebra of the system can also be used to 
establish (explicitly verifiable) criteria for dynamical realizability of unitary operators, 
as well as reachability /non-reachability of quantum states p2| , p^ . However, in general, 
controllability analysis does not provide a constructive procedure for realizing a desired 
control objective that has been shown to be theoretically achievable. Many techniques 
have been proposed to solve the challenging problem of finding suitable control fields. 
Most rely on some form of numerical optimization, model-based feedback or learning 



algorithms |2|, |5|, ||, gj, ||, ||, |3g, |3T|, |2 



The problem of constructive control using the Lie algebra structure of the system 
has also been addressed for certain systems such as coupled spin systems, where 
explicit procedures for generating arbitrary unitary operators have been developed 
|Tl|]. Constructive control techniques based on Lie group decompositions have also 



been developed to derive control schemes for explicit generation of unitary operators 



3^ , [35| , p6| |37| , |38| , and to achieve other control objectives such as maximization 
of observables and quantum state engineering in A^-level quantum systems with 
selectively addressable transitions. In similar techniques have been used to 



construct control sequences for a few target operators for some simple systems with 
degenerate transition frequencies and symmetrically coupled transitions. However, as 
regards constructive controllability for atomic systems with degenerate energy levels and 
transitions, further work is necessary. 

6.2. Controllability of open systems 

In this paper we systematically studied the degree of controllability of closed quantum 
control systems, i.e., quantum systems that only interact with a set of external control 
fields, whose interaction Hamiltonian is Hermitian. The question of controllability of 
open quantum systems, i.e., systems that also interact with the environment, is rather 
more complicated to answer. The time evolution operator of an open system depends 
critically on the interaction with the bath, represented by a non-Hermitian operator, 
and the dynamical evolution of the system is thus not unitary. Therefore, there are no 
conservation laws for energy, entropy and probability for such systems. 

This lack of unitary evolution and conservation laws for open quantum systems 
can be an advantage or an impediment. For instance, dissipation, e.g., by spontaneous 
emission, is essential for important applications such as optical pumping, quantum 



reservoir engineering with laser cooled trapped ions |41|, and laser cooling of internal 



molecular degrees of freedom |4^, |J]. Thus, one might say that dissipation enhances 
the controllability of these systems. On the other hand, there are applications such 
as quantum information processing, where dissipative effects such as decoherence are 
extremely deleterious. This suggests that the effect of dissipation on the degree of 
controllability of open quantum control systems should be assessed on a case-by-case 
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basis. 

7. Conclusion 

We have provided precise definitions for various degrees of controllability for closed 
quantum control systems as well as necessary and sufficient conditions for each, and 
applied the theoretical results to determine the degree of controllability of various atomic 
systems with degenerate energy levels and transition frequencies. Although the results 
of this controllability analysis based on the dynamical Lie algebra of the system often 
confirm our intuition about the controllability of the system, they can sometimes be 
surprising, such as the difference in the degree of controllability for transitions between 
states with the same degree of degeneracy, and those between states with different F- 
values. 

At present, our controllability analysis for atomic systems with degenerate energy 
levels and transition frequencies allows us only to determine the degree of controllability 
of the system, which can then be used to assess the feasibility of certain control objectives. 
In the future, it may be possible to use the structure of the Lie algebra of these systems 
to develop constructive techniques for the generation of unitary operators using simple 
control pulse sequences, as has been done for coupled spin systems and A'^-level systems 
with selectively addressable transitions. 
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Appendix A. Identifying the dynamical Lie algebra 

In order to decide whether the dynamical Lie algebra of a system with Hamiltonian (^) 
is of type sp(y), su{N) or u{N) we proceed as follows. We first compute the dimension 
of the Lie algebra generated by the operators iHm, i.e., the trace-zero parts of the iHm 
using the algorithm described in 0. If the dimension of the Lie algebra C is A^^ — 1 we 
check the trace of the operators Hm- If at least one of the operators has non-zero trace 
then the Lie algebra is isomorphic to u{N) and the system is completely controllable. If 
the dimension of C is A^^ — 1 and all Hm have zero trace then C = su{N) and the system 
is both density matrix and observable controllable. If is even and C has dimension 
A^(A^ + l)/2 then we check if there exists an operator J, which is unitarily equivalent 
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to the J defined in equation ( ]TT| ) such that 

{tHmfj + J{tHrn) = (A.l) 

holds for all m. If such a J exists then we conclude that the Lie algebra C is isomorphic 
to a subalgebra of sp^^) and as the dimension of the subalgebra equals the dimension 
of it follows that C ~ sp{y)- 

To determine if there exists a J such that ( [A.l| ) is satisfied we note that ( [A.l ) can 
be rewritten as system of linear equations 

£„J = 0, 0<m<M, 

where Cm is an A^^ x A^^ matrix determined by the iHm and J is an A^^ column vector. 
To determine the solutions J of the above equation we compute the null space of 

( ^ 

If it is non-empty then there is a J such that ([A.l|) is satisfied and we can compute its 



eigenvalues and compare them to the eigenvalues of the standard J given in (|TTD, which 
allows us to determine if this J is unitarily equivalent to the standard J. (Recall that 
two operators are unitarily equivalent if and only if the have the same eigenvalues.) 



Appendix B. Identifying forbidden operations 

To show that selective excitation of the m = sublevel population for a transition 
between two P3/2 states is indeed impossible in the context of this model, no matter what 
fields are employed, we note that this selective excitation would require implementing a 
unitary operator U that permutes the populations of states |3) and |4) (see Example 3a 
for notation) without disturbing the other states. Such a unitary operator would have 
to be of the form 
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modulo phase factors for the non-zero entries. However, it can easily be verified using 
the technique outlined in the previous section, that the Lie algebra generated by the 
system, i.e., by iHm with Hm as in equations ([T9|), (pO]) and (pi]), is isomorphic to sp{3) 
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Any unitary operator generated by this system must preserve this J, i.e., U^JU = J. It 
is easy to verify that the U above does not preserve J. Hence, it cannot be implemented. 
It is similarly easy to show that selective excitation of the m = — lorm = +l sublevel 
is not possible for this system. 



